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Abstract 



The results on polarized parton densities (PDFs) obtained using different meth- 
ods of QCD analysis of the present polarized DIS data are discussed. Their de- 
pendence on the method used in the analysis, accounting or not for the kinematic 
and dynamic 1/Q 2 corrections to spin structure function gi, is demonstrated. It is 
pointed out that the precise data in the preasymptotic region require a more careful 
matching of the QCD predictions to the data in this region in order to determine 
the polarized PDFs correctly. 



PACS numbers: 13.60.Hb, 12.38.-t, 14.20.Dh 



1 Introduction 



Our present knowledge of the spin structure of the nucleon comes mainly from polarized 
inclusive and semi-inclusive DIS experiments at SLAC, CERN, DESY and JLab, polarized 
proton-proton collisions at RHIC and polarized photoproduction experiments. One of the 
important and best studied aspects of this knowledge is the determination, in a QCD 
framework, of the longitudinal polarized parton densities and their first moments, which 
are related to the spins carried by the quarks and gluons in the nucleon. Different methods 
of analysis have been used in these studies. The aim of this paper is to discuss and clarify 
how the results on polarized PDFs depend on the method used in the QCD analysis. 

One of the peculiarities of polarized DIS is that more than a half of the present data 
are at moderate Q 2 and W 2 (Q 2 ~ 1 - 4 GeV 2 , 4 GeV 2 < W 2 < 10 GeV 2 ), or in the so- 
called preasymptotic region. So, in contrast to the unpolarized case, this region cannot be 
excluded from the analysis, and the role of the l/Q 2 terms (kinematic - 7 2 factor, target 
mass corrections, and dynamic - higher twist corrections to the spin structure function 
gi) in the determination of the polarized PDFs has to be investigated. This makes the 
QCD analysis of the data much more complicated and difficult than in the unpolarized 
case. 

2 QCD framework for inclusive polarized DIS 
2.1 Which data to chose for QCD fits 

The best manner to determine the polarized PDFs is to perform a QCD fit to the data 
on gi/Fi, which can be obtained if both the An and A±_ asymmetries are measured. 
In some cases only A\\ is measured. One can write (D is the depolarization factor and 
7 2 = 4MV/Q 2 ) 

^ = A l+V A 2 = (l + 7 2 )|- + ( V - l)A 2 , (1) 

from which one sees that the quantity Au/D(l + 7 2 ) is a good approximation of g\jF\ 
because the second term in the second relation of (1) can be neglected in the preasymptotic 
region too - the asymmetry A 2 is bounded and in fact small, and multiplied in addition 
by a small kinematic factor (77 7) . 

The data on the photon-nucleon asymmetry A\ are not suitable for the determination 
of PDFs because the structure function g 2 is not well known in QCD and the approxima- 
tion 

= g 1 /F 1 - 1 2 g 2 /F l » (g 1 /F 1 ) theor (2) 
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used by some of the groups is not reasonable in the preasymptotic region because 7 2 cannot 
be neglected. 

Bearing in mind the remarks above, let us discuss in more detail how to confront 
correctly the theoretical predictions to the available polarized inclusive DIS data: 

i) First of all, one should include in the QCD fit of the world data all g\/F\ data 
available. These are the CLAS(p, d), JLab/Hall A(n), SLAC E143(p, d) and E155(p, d) 
data [l]. 1 

ii) For the rest of experiments: EMC(p), SMC(p, d) [2] and COMPASS(d) [3] at 
CERN, HERMES (p, d) [4] at DESY and E142(n), E154(n) at SLAC [5], only data on 
A 1 are presented. In the experiments at CERN and DESY only the asymmetry A\\ is 
measured. However, for different reasons the approximations A\\/D ~ A\ ps gi/Fi for 
CERN data and A\\/[D(1 - 777)] w A 1 w g 1 /F 1 for HERMES data are good ones. For the 
experiments at CERN, the 7 2 factor is very small and the term / y 2 g2/F 1 can be neglected, 
while for the HERMES data the approximation g 2 = is used and the effect of the non- 
zero value of g 2 is included in the systematic uncertainty of Ai. 2 In the SLAC E142(n) 
and E154(n) experiments both A\\ and A± have been measured and gi/F-y data could 
have been extracted, but the collaborations present only data on A\. Bearing in mind 
the kinematic region (E154) and the precision (E142) of these data, the approximation 
Ai/ (I +7 2 ) ps gi/Fi in Eq. (1) for them is reasonable. 

To summarize, in the pure DIS region A 1 ps g i /F 1 and it does not matter which 
data are be used in the QCD analysis. This is not the case when precise data in the 
preasymtotic region have to be used too. In that case one has to confront the QCD 
predictions to the data more carefully in order to extract the polarized PDFs correctly. 

2.2 Methods of QCD analysis 

In QCD, one can split g\ and Fi into leading (LT) and dynamical higher twist (HT) pieces 

<?i = (<?i)lt,tmc + (a , i)ht, F\ = (Fi)lt,tmc + (Fi)b.t- (3) 
In the LT pieces in Eq. (3) the calculable target mass corrections (TMC) are included 

gi(x, <5 2 )lt,tmc = gi(x, Q 2 )lt + gi(x, Q 2 )tmc, (4) 

1 Note that excepting the E155 Collaboration, the other Collaborations present data on A\ too. The 
corresponding values of A\ and g\jF\ at the same (x, Q 2 ) are different, which means that the term 
r ) 2 g2/F\ cannot be really neglected in the preasymptotic region and a fit to A\ data instead of (gi/F\), 
approximating (Ai) theor with (gi/Fi) theor is not correct. 

2 Note that for the final inclusive HERMES A\ data [6] the approximation g 2 = is not used and the 
relation A\ « gi/F\ does not hold. 



2 



Fl(x,Q 2 ) LT , T MC = F 1 (x,Q 2 ) LT + F 1 (x,Q 2 ) T MC- 

They are inverse powers of Q 2 kinematic corrections, which, however, effectively belong 
to the LT part of Q\. Then, approximately 

9l ^ (gyiST ^ _|_ (ffpTMC+HT _ (-^i)tmc+ht j ^ 
Fi (-Fi)lt Q?i)lt (-Pi)lt 

Note that the LT pieces (<7i)lt an d (-F)lt are expressed in terms of the polarized and 
unpolarized PDFs, respectively. In what follows only the first terms in the TMC and HT 
expansions will be considered 

gi(x,Q 2 ) TMC = M 2 /Q 2 g^\x 1 Q 2 ) TMC + 0(M A /Q i ) 1 (6) 
Fi(z,Q 2 )tmc = M 2 /Q 2 Fi 1 \x,Q 2 ) TMC + 0(M 4 /Q% 

9i(x,Q 2 ) RT = h^(x)/Q 2 + 0(A 4 /Q 4 ), 2xF 1 (x,Q 2 ) RT = h 2xF ^x)/Q 2 + 0(A 4 /Q 4 ). (7) 

The first terms of the HT pieces in Eq. (7) are shown in Fig. 1 and as seen, they are 
definitely different from zero and cannot be neglected in the preasymptotic region. 



rS - ' 0.1 

> 

o 

-x 



-0.2 



-0.3 



proton 



T? * * + 



a * m 2 



■ g^LSS'OB) 
o F 2 (Al'02) 
• 2xF 1 (Al'02) 



o.o 



0.2 



0.4 



0.6 0.8 
X 



Figure 1. HT corrections to g\ [10], i*2 and 2x.Fi [11] structure functions. 



There are essentially two methods to fit the data - taking or NOT taking into account 
the HT corrections to g\. According to the first [7], the data on g\/F\ have been fitted 
including the contribution of the first term h(x)/Q 2 in (#i)ht and using the experimental 
data for the unpolarized structure function Fi 



9i(x,Q 2 
Fi{x,Q 2 



gi(x,Q 2 ) LT ,TMc + h^(x)/Q 2 



exp 



f^q 2 : 



(Method I). 
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According to the second approach [8] only the LT terms for g\ and F\ in (3) have been 
used in the fit to the gi/Fi data 



9i(x,C?)] ^ g(;^2 LT (Method II). (9) 
_ exp F 1 (x,Q 2 ) LT 



Fi{x,Q 



2\ 



It is obvious that the two methods are equivalent in the pure DIS region where HT 
can be ignored. To be equivalent in the preasymptotic region requires a cancellation be- 
tween the ratios (#i)tmc+ht/(#i)lt and (Fi)tmc+ht/(-Pi)lt in (5)- Then (#i)lt obtained 
from the best fit to the data will coincide within the errors independently of the method 
which has been used. In Fig. 2 these ratios based on our results on target mass [9] 
and higher twist [10] corrections to g±, and the results on the unpolarized structure func- 
tion Fi are presented. Note that for the neutron target, ((?i)tmc+ht is compared with 
(Pi)lt ^ 1 (j^)lt HT because of a node-type behaviour of (<7i)lt- Also, LT means the NLO 
QCD approximation for both, gi and F 1 . As seen from Fig. 2, (TMC+HT) corrections 
to gi and F\ in the ratio g\/F\ do not cancel and ignoring them using the second method 
is incorrect and will impact on the determination of the polarized PDFs. 3 




Figure 2. Comparison of the ratios of (TMC+HT)/LT for g\ and F\ structure functions for 
proton and neutron targets (see the text). 

Modifications of the second method of analysis in which F\ is treated in different way 
are also presented in the literature: 

3 Note that this result differs from our previous observation that the ratios (<7i)ht/(<7i)lt and 
(-Fi)ht/(-Fi)lt for a proton target approximately cancel for x > 0.15 at Q 2 = 2.5 GeV 2 [12]. How- 
ever, in a more precise analysis one should account for the TMC corrections in Eq. (5). Note that 
previously, for the calculation of the ratio (i 7 i)HT/(^ 7 i)LT the results of [11] were used, where a QCD 
analysis of the world unpolarized DIS data was performed using the cut Q 2 > 2.5 GeV 2 . This cut ex- 
cludes a lot of data in the preasymptotic region which influence the determination of the HT corrections 
to F\. Fig. 2 is based on a new analysis of the ratio (FiJtmc+htA-FiJlt [13]. 
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i) Blumlein, Bottcher [14] and COMPASS [3], where instead of (Fi) LT in (9), (Fi) 
has been used in the fit to the data 

~ 9i{x,Q 2 ) 
_Fi(x,Q 2 ) 

ii) AAC Collaboration [15], where Fi is expressed in terms of F 2 and R, and for them: 
(-^2)lt and R exp have been used, respectively 

Ai(x, Q 2 )exp ~ 

As mentioned above, the approximation for Ai in (11) is not correct for most of the data 
sets used in the fit. 

Note that when the second method or its modifications are used, the HT effects of gi 
are apparently absorbed into the extracted PDFs, which thus differ from those determined 
in the presence of HT (for more details see the discussion below), but, of course, not all 
the data can be fitted satisfactorily. On other hand, the extracted PDFs in the framework 
of the second method and corresponding to fits (9 - 11) should all be different due to the 
HT corrections to F\ (see Fig. 1) which are included in Eq. (10) but not in Eq. (9), and 
only partly in Eq. (11). So, using the different denominators in (9 - 11) one will obtain 
different values for the free parameters associated with the input polarized PDFs after 
fitting the data. 

3 Polarized Parton Densities 

We will discuss in this section in more detail how the results on polarized PDFs depend on 
the method used for their determination. To illustrate this dependence we will compare 
the NLO LSS'06 set of polarized parton densities [10] determined by Method I with those 
obtained by COMPASS [3], DSSV [16] and AAC Collaboration [17] using Method II or 
its modifications. 

3.1 Comparison between LSS'06 and COMPASS PDFs 

To obtain the LSS'06 PDFs we used Method I (Eq. 8) in the QCD analysis of the the 
world data on polarized inclusive DIS ([l]-[5]), i.e. the HT corrections to g\ were taken into 
account. For the LT term we have used the NLO QCD approximation in the MS renor- 
malization scheme. In their analysis COMPASS has used Eq. (10), but the CLAS data 



t rn 



gi(x,Q 2 ) LT 

F\ (x, Q 2 ^)exp 



(10) 



9i(x,Q 2 



;n) 
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from [1] were not included in their fit. 4 In both the analyses the experimental data for the 
unpolarized structure function F\ were used. Thus what is fitted by (gf^ur^COMPASS) 
is significantly different from what is fitted by our (gf^LT^LSS), i.e. 

g?(x, Q 2 )lt(COMPASS) = g?(x, Q 2 )lt,tmc(LSS) + h N (x)/Q 2 = ( f) t L o s t s (N = p, n, d). 

(12) 

Note that for (<?i)lt COMPASS has also used the NLO QCD approximation in the MS 
scheme. As a result, (gi)^ 3 an d (<?i)lt(COMPASS) obtained from the fit are almost 
identical, but the LT terms of g\ corresponding to LSS and COMPASS fits are different 
for x < 0.1 where the HT corrections to gf cannot be neglected. This fact is illustrated 
in Fig. 3 for gf where the COMPASS data are also presented. We have found that the 




Figure 3. Comparison between (gi)tot(LSS) and {g\)LT{COMPASS) obtained from the fit 
with the COMPASS data at measured x and Q 2 . Error bars represent the total (statistical and 
systematic) errors. The (gi)LT{LSS) curve is also shown. 

HT contribution to (<?i)tot, h d (x)/Q 2 , is positive and large, up to 40% of the magnitude of 
(gf)vr in the small x region, where Q 2 is small (Q 2 ~ 1 — 3 GeV 2 ). As a consequence, the 
HT effects are effectively absorbed in the COMPASS PDFs. A crucial point is that the 
COMPASS analysis does not include CLAS data, which are entirely in the pre-asymptotic 
region, and for which the HT effects are essential. In Fig. 4 the COMPASS PDFs 
corresponding to a positive solution for AG are compared with those obtained by LSS'06. 
As seen from Fig. 4, except for (Aw + Au) the other PDFs differ, especially those of the 
strange quarks and gluons. 

4 Using this method one cannot achieve an acceptable value of x 2 for the CLAS data which are entirely 
in the preasymptotic region [18]. 
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Figure 4. Comparison between NLO(MS) LSS'06 polarized PDFs and those obtained by 
COMPASS. 

3.2 Comparison between LSS'06 and DSSV PDFs 

Recently the DSSV group has presented results on NLO(MS) polarized PDFs [16] obtained 
from the first global analysis of polarized DIS, SIDIS and RHIC polarized pp scattering 
data. Due to the SIDIS data a flavor decomposition of the polarized sea is achieved. 
For the fit to the inclusive DIS data the second method (9) was used, i.e., a NLO QCD 
approximation for ((?i)lt and (-Fi) L t in the ratio g\fF\. The unpolarized structure function 
Fx(x, Q 2 )lt was calculated using the NLO MRST'02 parton densities [19]. The difference 
between Fi(x, Q 2 )nlo and the phenomenological parametrization of the data, F\(x, Q 2 ) exp , 
used in our analysis [10] is illustrated in Fig. 5. It is a measure of the size of the TM and 
HT corrections Fi(x, Q 2 )tmc+ht to F\ which cannot be ignored in the pre-asymptotic 
region. Note that in the MRST fit to the unpolarized data the preasymptotic region 
was excluded precisely in order to eliminate the TM and HT corrections. That is why 
-Fi(MRST)nlo differ from (-Fi) ex p in the preasymptotic region where the TM and HT 
corrections to F\ cannot be ignored. 

It is important to mention also that in the preasymptotic region for large x and lower 
Q 2 (x > 0.40, 0.47 for JLab and SLAC/E143 data, respectively), the kinematic factor 
7 2 = 4M 2 x 2 /Q 2 is larger than R(x,Q 2 ). Then it follows from the relation between F\ 
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Figure 5. Comparison between Ff (MRST)nlo an d (Ff) exp unpolarized structure functions at 
Q 2 = 1.5 GeV 2 (left) and Q 2 = 2.5 GeV 2 (right). 



and F 2 



2xF 1 (x,Q 2 ) =F 2 (x,Q 2 



1+7 2 



l + R(x,Q 2 ) 

that in this region 2xFi > F 2 and, as a consequence, 

F L (x, Q 2 ) = F 2 (x, Q 2 ) - 2xF 1 (x, Q 2 ) 



(13) 



(14) 



is negative, in contrast to what follows from pQCD, i.e. that (F L ) LT should be always 
positive. So, F L could become negative in the preasymptotic region due to HT corrections. 
We consider it is important to test this observation by fitting the data on unpolarized 
structure functions F 2 and Fl which will become available at JLab in the near future. 

The main features of the results of the fits obtained by LSS (Method I) and DSSV 
(Method II) are illustrated in Fig. 6 for a proton target. As expected, the curves cor- 
responding to the ratios <7*°*(LSS)/(Fi) cxp and pi(DSSV)nlo/Fi(MRST)nlo practically 
coincide although different expressions were used for g% and Fi in the fit (see the left 
panel of Fig. 6; the difference between them for x > 0.2 will be discussed later). In the 
right panel of Fig. 6 the LSS and DSSV LT(NLO) pieces of g\ are compared for a proton 
target. Surprisingly they coincide for x > 0.1 although the HT corrections, taken into 
account in LSS'06 and ignored in the DSSV analysis, do NOT cancel in the ratio gi/F\ 
in this region, as has already been discussed above. The understanding of this puzzle is 
connected with the fact that in the DSSV fit to all available gi/F% data a factor (1 + 7 2 ) 
was introduced on the RHS of Eq. (9) 



gi{x,Q< 



Fi(x,Q 2 ) 



9i(x, Q 2 )vt 



J exp 



;i + 7 2 )F 1 (a;,Q 2 ) L T 



(15) 



(Note that for the fit to the A\ data Eq. (9) was used. 
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Figure 6. Comparison between: the ratios ^ ot (LSS)/(Fi) exp and #i(DSSV) N lo/Fi(MRST) nlo 
(left); <?i(LSS)nlo an d 5'i(DSSV)nlo (right). The LSS results correspond to the node-type 
solution for xAG(x,Q 2 ) . 

There is no rational explanation for such a correction. The authors point out [20] that 
it is impossible to achieve a good description of the gi/F\ data, especially of the CLAS 
ones, without this correction (see Fig. 7). As seen from Fig. 7, the theoretical curves 
lie systematically above the data which are badly fitted without introducing the (1 +7 2 ) 
factor. 




Figure 7. CLAS g\/F\ data compared to the theoretical DSSV curves accounted or not for the 
(1 + 7 2 ) factor. 

It turns out empirically that the 1/(1 + r y 2 ) factor accidentally more or less accounts 
for the TM and HT corrections to g\ and F\ in the ratio g\j F\. The relation 

1 (gpTMC+HT (Jl)TMC+HT 1 

(0l)lT " (fl)LT ^(l + 7 2 ) 1 ] 



9 



is satisfied with an accuracy between 4% and 18% for the CLAS proton data (x > 0.1 and 
Q 2 between 1 and 4 GeV 2 ). That is the reason why the LT(NLO) pieces of gf obtained by 
LSS and DSSV are in a good agreement for x > 0.1 (see the right panel of fig. 6). Also, 
why the curve in Fig 6 (left) corresponding to (7i(DSSV)nlo/Fi(MRST)nlo lies above the 
one of g\ ot (LSS) / (Fi) exp . Including the (I+7 2 ) factor in (15) would make the curves almost 
identical. It is important to mention that introducing the (1 + 7 2 ) factor does not help 
at x < 0.2 because it is small in this region and cannot mimic the difference between TM 
and HT corrections to g\ and F\ (LHS of Eq. (16)) for proton as well as for neutron target 
(see Fig. 2). That is why the LT(NLO) pieces of g\ obtained by LSS and DSSV groups 
differ in this region - the smaller x is, the greater is the difference (see Fig. 6 (right)). 
To summarize: It is impossible to describe the precise data in the preasymptotic region 
like the CLAS data using the second method. Its empirical modification by introducing 
the (1 + 7 2 ) factor accounts approximately for the TM and HT effects, but only in the x 
region: x > 0.1, 0.2 for proton and neutron targets, respectively. 

The NLO(MS) PDFs determined from LSS'06 and DSSV analyses are compared in 
Fig. 8. The AAC'08 PDFs [17] obtained from a combined NLO QCD analysis of inclusive 
DIS and RHIC 7r°-production data are also presented. Note that for the fit to DIS data 




0.03 

0.02 
0.01 
0.00 
■0.01 
-0.02 

■0.03 



I 

■ xAs Q 2 = 2.5 GeV 2 


l 

\ \ 
\ 1 




\ \XS(X) 




\ \ 












AT 




/ / 1 

y s 1 


MRST'02 


1 1 


GRV'98 


1 1 



0.01 



0.1 




Figure 8. Comparison between LSS'06, DSSV and AAC'08 NLO PDFs in (MS) scheme. 

AAC have used the modification (11) of the second method and that the RHIC data 
impact mainly on AG. Note also that all these analyses include the precise CLAS data in 
the preasyptotic region. The results are presented for the sums (Au + Au) and (Ad + Ad) 
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because they can only be separated using SIDIS data (the DSSV analysis). Although the 
first moments obtained for the PDFs are almost identical, the polarized quark densities 
themselves are different, especially As(x) (in all the analyses As(x) = As(x) is assumed). 

Let us discuss the impact of HT effects on (Au + Au) and (Ad + Ad) parton densities 
which should be well determined from the inclusive DIS data. (Au + Au) extracted by 
LSS and DSSV are well consistent. As was discussed above the HT effects for the proton 
target are effectively accounted for in the DSSV analysis for x > 0.1 by the introduction 
of the (I + 7 2 ) factor. However, this factor cannot account for the HT effects for a neutron 
target at x < 0.2 (see Fig. 2) and the impact of higher twist on (Ad + Ad) determined by 
DSSV is demonstrated in Fig. 8. The positive HT effects are absorbed into (Ad+ Ad) D ssv 
and it is thus less negative in this region. The influence of HT effects at small x (not 
accounted for by the DSSV group) on both parton densities is not sizable because of two 
reasons: first, their values are small and second, the data in this region are not precise 
enough to indicate the impact of higher twist on their values. The impact of HT effects on 
both (Am + Am)aac and (Ad + Ad) AAC is larger because the AAC Collaboration has not 
taken them into account at all, and in addition, the incorrect approximation A\ g\jF\ 
for some of the data in the preasymptotic region has been used. The difference between 
the strange sea densities As(x, Q 2 )lss an d As(x, Q 2 )aac f° r a; > 0.1 is due to the different 
positivity conditions which have been used by the two groups. Note also that the positivity 
condition \AG(x, Q 2 )\ < G(x, Q 2 ) is not satisfied for the polarized gluon density obtained 
by AAC, which suggests it is not physical. 

In contrast to a negative As(x, Q 2 ) obtained in all analyses of inclusive DIS data, 
the DSSV global analysis yields a changing in sign As(x, Q 2 ): positive for x > 0.03 
and negative for small x. Its first moment is negative (practically fixed by the SU(3) 
symmetric value of a 8 ) and almost identical with that obtained in the inclusive DIS 
analyses. It was shown [21] that the determination of As(x) from SIDIS strongly depends 
on the fragmentation functions (FFs) and the new FFs [22] are crucially responsible for 
the unexpected behavior of As(x). So, obtaining a final and unequivical result for As(x) 
remains a challenge for further research on the internal spin structure of the nucleon. 

4 Summary 

The fact that more than a half of the present polarized DIS data are in the preasymptotic 
region makes the QCD analysis of the data more complex and difficult. In contrast to 
the unpolarized case, the 1/Q 2 terms (kinematic - r j 2 factor, target mass corrections, and 
dynamic - higher twist corrections to the spin structure function gi) cannot be ignored, 
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and their role in determining the polarized PDFs is important. Sets of polarized PDFs 
extracted from the data using different methods of QCD analysis, accounting or not 
accounting for the kinematic and dynamic 1/Q 2 corrections, are considered. The impact 
of higher twist effects on the determination of the parton densities is demonstrated. It is 
pointed out that the very accurate DIS data in the preasymptotic region require a more 
careful matching of QCD to the data in order to extract the polarized PDFs correctly. 
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